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SUMMARY 


Flight tests and a mathematical analysis were made to demonstrate 
and confirm a type of suhsonlc flutter involving rigid-body motions and 
ving deformations. For the configuration considered, the period of the 
oscillation was approximately 100 chords per cycle which is well within 
the range of period found in dynamic- stability work on rigid aircraft 
with free controls. A mathematical analysis based on two-dimensional 
incon^jres Bible flow provided a conservative prediction of the airspeed 
at which the low-frequency flutter occiorred. It was fo\md that wing 
bending stiffness is the inq)ortant pareuaeter for preventing such flutter. 


INTRODUCTION 


Interaction of deformations of an aircraft structure with the passing 
alrstream can lead to the dynamic instability known as flutter. For 
bending-torsion wing flutter, the frequency of oscillation is fairly high 
and usually approaches the natural torsional frequency of the wing in 
still air. Such an oscillation may be contrasted with ordinary dynamic- 
stability phenomena involving rigid-body modes which lead to much lower 
frequencies that are usually controllable. 

The fact that the calculated flutter speed may be modified by the 
addition of free-body modes has been recognized for many years . For 
exanqjle, about 20 years ago it was found analytically (ref. l) that body 
mobility had a slight favorable effect on the calculated flutter speed 
of a parti (nilar configuration typical of that day. The problem has from 
time to time been reconsidered in both American and British literature 

^Supersedes the recently declassified NACA RM L50I29, "description 
and Analysis of a Rocket- Vehicle Ejgieriment on Flutter Involving Wing 
Deformation and Body Motions" by H. J. Cunningham and R. R. Lundstrom, 

1950. 
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and the necessity for determining any potential detrimental effect of 
special configurations, sweptback wings, higher speeds, and higher alti- 
tudes has lately become more insistent. A recent paper by Broadbent 
(ref. 2) discusses the necessity for including free-body modes in the 
study of sweptback wings. 

Controlled experimentation involving free-body modes is highly desir- 
able although difficult even at low speeds in a wind tunnel. In refer- 
ence 5 Lambourne gives some experimental results and describes various 
difficulties encountered, principally the difficulty of supporting the 
model so that act\ial flight behavior is sufficiently well simulated. 

The National Advisory Committee for Aeronautics has been engaged for 
several years in making flutter experiments in the transonic range by 
the \ise of modern telemetering techniques with rocket-powered vehicles. 

Some of the results of these experiments are given in several NACA papers 
Including references 4 and 5* Such flutter research has been concerned 
almost entirely with wing bending- torsion flutter. Recently, however, in 
connection with some of these experiments, flutter failures were obtained 
which definitely Involved much lower frequencies than those obtained in 
ordinary bending- tors ion flutter . A preliminary account of an unexpected 
low-frequency failure is given in reference h, and two other low-frequency 
failures have been observed by W. T. Lauten, Jr., and J. M. Teitelbam 
during rocket-propelled model tests at the Langley Pilotless Aircreift 
Research Station at Wallops Island, Va. 

It was decided to repeat the e:q)erlment of the D model of reference 4 
with more extensive instrumentation in order to obtain information specific 
to this type of flutter. The present paper includes a description of and 
results from model D of reference as well as similar and more conprehen- 
sive material for the repeat ejqjeriment on the model which has been desig- 
nated model E. Resiilts of an aneOysls developed on the basis of certain 
simplifying assuoptlons are given, and a ccxnparlson is made between analyt- 
ical and experimental results. 


SYMBOLS 


a 


nondimenslonal position of wing elastic axis, - 1 


100 


a + x^ nondimenslonal position of wing- section center of gravity. 



100 



geometric aspect ratio of one exposed wing panel, for rectangular 
plan-form wings, Length/Chord 
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b wing semi chord, ft 

El bending rigidity, Ib-in.^ 

B nondimensional distance (in ving semi chords) from wing midchord 

line to missile center of gravity meas\xred parallel to missile 
longitudinal axis, positive when center of gravity is behind 
wing midchord line 

fh^ first-bending natural frequency, cps 

fhg second-bending natural frequency, cps 

ft first-torsion natural frequency, cps 

fa uncoupled first- torsion frequency relative to elastic axis, 

1/2 

, cps 

structural d a m p ing coefficient in bending 

£!^ structural damping coefficient in torsion 

g mathematical quantity having qualities of structiiral daaqplng 

coefficient 

GJ torsional rigidity, Ib-ln.^ 

lyy mom ent of inertia of missile about center of gravity, slug-ft^ 

k reduced frequency parameter, cob/v; jt/k is period of oscillation 

in chords per cycle 

Z length of wing along leeuiing edge outboard of body, in. 

m mass of wing per unit length, slugs /ft 

M Mach number 

theoretical Mach number at which sonic velocity is first attained 
over section of wing taken perpendicular to leading edge at 
zero lift 


1 - 






k 
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■•a 


nondlmensional radius of gyration of wing section about elastic 


V 

bcuh 

^cg 

Xea 


K ( stnd • ) 
P 


ajcis 


i£/ 


mb'^ where 3^ is the mass moment of inertia of 
wing about its elastic axis per unit length in slug-ft^ft 

time after rocket laimching, sec 
speed, ft/sec 

nondlmensional flutter-speed coefficient 

distance of center of gravity of wing section behind leading 
edge, percent chord 

distance of wing elastic sals behind leading edge, percent 
chord 

mass ratio, «pb^ m 

value of K when p is standard sea- level density 
air density, slugs/cu ft 


CO 

angular 

“'h 

angular 


angular 

2rtfa,, 


flutter 


flutter parameter, (l + 


EXPERIMENTAL INVESTIGATION 
Apparatus and Methods 


The test vehicle was essentially a tailless configuration having the 
test wings as the only stabilizing surfaces in pitch. The model was 
powered by a modified Aerojet 12AS-1000 D rocket motor capable of carrying 
it to greater- than- sonic speed with an acceleration of about times grav- 
ity. A sketch of model E is shown in figure 1. Center of gravity, weight, 
and moment of inertia in pitch varied with flight time because of the fuel 
consumed. These parameters plotted against flight time for model E are 
shown in figure 2. The laimchlng was made at the Langley Pilotless Aircraft 
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Research Station at Wallops Island, Va., from a near- zero- lei^h launcher, 
as descrihed in reference 5* Photographs of model E in the launching 
position are shown in figure 5- The test wings were made of laminated 
white pine with an inlay of 0.032-inch 2i*S-T aluminum alloy to duplicate 
as closely as possible the test wings of model D of reference 4. Strain 
gages were mounted on the elastic axes of the wings 4 inches from the wing 
root, as shown in figure 1. 

A six- channel telemeter was installed in the model with instrumen- 
tation to give continuous readings of the strains, left-wing bending, 
right-wing bending, and right-wing torsion, at the coiresponding gage 
locations and also to give continuous readings of longitudinal and normal 
accelerations of the model center of gravity and angle of attack of the 
missile. Speed of the model was obtained by integration of longitudinal 
acceleration, and the altitude was obtained from a pxilse-type tracking 
radar \mlt. Atmospheric conditions prevailing at the time of flight were 
obtained from a radiosonde. 

The accuracy of the quantities meastired is believed to be within the 
following limits; longitudinal and normal accelerations, 2 percent; vel- 
ocity, 4 percent; wing bending and twisting moments, 20 percent; and phase 
angle between any two quantities, 15 °. 


Experimental Results 

The physical characteristics of the wings as determined by preflight 
groimd tests are listed in table I. Data for model E are presented in 
figures 4 to 7 as functions of the flight time. Included in these figures 
are the available data for model D from reference 4 for comparison. Fig- 
ure 4 includes flight velocities; figure 5 shows Mach number and air density; 
figure 6 shows longitudinal acceleration. Figure 7 presents the following 
data; bending moments of both wings (positive, wing bent down) at the 
strain-gage locations, twisting moment of the right wing about its elastic 
axis (positive, leading edge up) at the strain-gage location, and normal 
acceleration of the missile center of gravity (positive, up) . It was 
desired to obtain the angle of attack of model E directly as a function 
of time, but the angle- of -attack indicator was inoperative dioring flight. 

As a result, the angle of attack was determined by use of the normal accel- 
eration and an assumed lift-ciirve slope of 0.08 per degree as obtained 
from lifting- surface theory. 

Comparison of figure 7 for model E and figure 9 of reference 4 for 
model D indicates that both flights had similar behavior. At a speed of 
about 400 feet per second, there began pitching oscillations with slowly 
increasing amplitudes which continued to a speed of about 620 feet per 
second for model D and 500 feet per second for model E. Over this inter-val 
the primary frequencies of oscillation increased from 4 to 7 cycles per 
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second. for model E and from 5 to 8 cycles per second for model D, and 
there was a considerable component of the wing first-bending natural fre- 
quencies present on the wing bending strain-gage traces. An interruption 
or danping of the oscillations then occurred and was followed by final 
divergent oscillations with frequencies of about 9*0 cycles per second 
for model D and 8.5 cycles per second for model E which resvilted in fail- 
vtres of the wing after about 7 or 8 cycles. The final oscillations 
were of a much more regular character and approached a ^vergent sine 
wave. The speeds at the beginning of the final divergent oscillations, 

670 feet per second for model D and ^dO feet per second for model E, are 
considered to be the experimental flutter speeds beca^lse of the regularity 
of the subsequent oscillations. As indicated by the time histories, the 
periods of the low-frequency oscillations for both models, from the begin- 
ning of a detectable oscillation until wing failure, were about 100 chords 
per cycle. 

The primary interest in flutter research centers on the speed end 
frequency of oscillation. A secondary interest exists in the deflection 
amplitudes and in the phase relationships of the varloiis motions j the 
latter are usually not experimentally determinable with good accuracy. 

From analysis of the oscillograph records of model E taJcen during flutter 
over the interval t = 5*^0 seconds to 5*65 seconds, the amplitudes of 
wing deformation and body motions relative to a wlng-tlp-bending amplitude 
of 1 inch were found to average about O.OO65 radian for wing- tip torsion, 
0.020 radian for pitching, and 0.21 inch for vertical translation. The 
wing-tlp-bending and wing- tip-torsion amplitudes were determined through 
the assumption of a constant lift distribution from wing tip to wing tip, 
whereas the amplitude of vertical translation was obtained by integrating 
the normal acceleration twice. The body-pitching anplitudes were determined 
by using the normal acceleration and an assumed lift-curve slope of 
0.08 per degree. The overall accxuracy of the experimentally determined 
amplitudes is thought to be within ±50 percent^ thus, these values are 
at least of the proper order of magnitude. 

It can be seen from figure 7 that the phase relationship is as fol- 
lows within the limit ±15° of experimental accuracy: At 8ui instant when 

the wing is bent up a maximum amount, it is also twisted (leading edge 
up) a maximum amount and, as Interpreted from the normal acceleration, 
the missile is pitched nose up and translated down a max1 imm amount. 


ANALYTICAL INVESTIGATION 
Bases of Analysis 


Only an outline of the analytical treatment vised to obtain the results 
is given here. It is assumed that the theory of linear superposition for 
small disturbances holds. A simplified configuration which has four degrees 


NACA TN 5511 


7 


of freedom in a symmetric oscillation is treated. There are two rigid- 
body degrees of freedom; pitching about the center of gravity, and ver- 
tical translation^ and two wing-deformation degrees of freedom which are 
approximated by the first-uncoupled-bending and first-uncoupled- torsion 
modal shapes of an ideal uniform cantilever beam. For convenience, the 
aerodynamic coefficients used for a wing section are those of incon 5 )res- 
sible, two-dimensional potential flow and yield forces which are propor- 
tional to the displacement and motion of that section (strip analysis) . 
Spanwise distributions of air forces due to symmetric motions of a rigid 
body are thus constant in magnitude and phasing from wing tip to wing tip, 
including the body intercept. In view of the Mach nximber of the experi- 
mental flutter, -vdiich is well below the critical Mach n\miber Mqj. of the 

airfoil section, aerodynamic coefficients for incoD5)reBsible flow eire 
considered satisfactory. Coefficients for compressible flow could be 
substituted, however. 

The configuration has no other horizontal airfoil surface and no 
other aerodynamic forces Eire assumed on the body. Each half -wing is 
uniform and the two half-wings, ri^t and left, of each missile are 
treated as being identical. The effect of gravity on an oscillation 
having the frequency of the observed oscillation is considered negligible, 
as are the effects of rocket thrust. 


Equations of Equilibrium 

In order to obtain the theoretical conditions for flutter, Lagrange's 
equation is employed, as, for example, in reference 6, to derive four 
equations of equilibrium. With the use of the specified condition of 
hEirmonic motion, the four equations may be written in matrix form as 
follows : 


ail a ]_2 ai 3 a^^ 


*1 



a2i a22 »25 a24 


X2 


c; 

aji aj2 a54 


X5 


ftpx^ 

_aiil ®45 ®’44 


3- 


OTx4 


2 

In this equation the qviantity (l + ig) where g has the prop- 

erties of a structijral damping coefficient. (Structural damping force 
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is proportional to amplitude and in phase with velocity.) The quantities 
xi, X 2 , xj; and represent the nondimens ional amplitudes of oscil- 

lation in wing bending, wing torsion, body pitching^ and vertical trans- 
lation, respectively. On the basis of the foregoing assumptions, for 
the configuration treated. 


an = 
ai2 = 
ai5 = 
ai4 = 

s-21 = 
&22 = 
^25 = 


K ■ 

1.5558(^ - 
1 .; 662 ( -^ a . Xa _ 

i.5662an 



1.2752 



+ 


(-S + 

K 



a^i^ = 1.8784a2i 
where 


'51 


= 0.55177 (~^ + a + xa _ 


a^2 = 0.5595 


^ (-S + a)xg, _ 


Ada 


'55 


- A(3^p 


= '-^dh 


a4i = 0.55177an 


a^2 = 0.59795aj_2 


%5 ~^cp 

^44 = - -^ch 


^ch “ ” "k" ^ k 


a H- 2F _ 


a>iiG + i p 

L 20 + 

(- - 

k2 

\2 

) k E 

c k2 

\2 jkj 


Aah = Aah(k,a) = - | - ^| + a^Ach 
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Aqp - 

Adh = Aah(k,s) 
Adp = Aaa(k, s) 


Aap(^>®')®) — Ann. “ (s - 8’)A q'|- | 
A^(k,a,s) = Ago, - (s - a)A^ 


and 


_ 2 

^p _ Mass moment of missile in pitch about its center of gravity 
K (Span)xpb^ 


1 _ Total mass of missile 
(Span)jfpb2 


The quantity a is the nondimens lonal distance (in wing semichords) 
from the wing midchord line to the wing elastic axis, positive when the 
elastic axis is behind the midchord. The quantity s is the nondimen- 
sional distance (in wing semichords) from the wing midchord line to the 
center of gravity of the missile, positive when the center of gravity is 
behind the midchord. The quantities F and G are the real and imagi- 
nary parts, respectively, of the complex function C = C(k) = F(k) + iG(k) 

developed in reference 7* The reduced frequency parameter k equals 

It is to be recognized that equation (l) represents a characteristic- 
value problem of aeroelastic harmonic motion. The depend on the 

reduced frequency k, whereas ft contains the frequency cj. Combinations 
of o) and k which result in a specified value of g are the character- 
istic values (eigenvalues) of the system. Knowledge of k and o) leads 
to a flutter speed v. The trivial solution xq = X 2 = xj = x4 = 0 is, of 
coxorse, not sought . 
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The quantities p and T are proportional to any mechanical 
restoring force due to body motions in the hody-pltching and vertical- 
translation degrees of freedom, respectively. Such restoring forces could 
exist, for example, if the model were moun ted on springs in a wind t^mnel . 
When the wing-body combination is of the freely flying type, however, there 
is no mec h a ni cal restoring force; that is, the natviral (in vacuo) frequen- 
cies of the two rigid-body degrees of freedom are zero. Hence, the 
characteristic- value solutions for equation (l) are sought at the 
limit as p and T approach zero. 


DISCUSSION AND COMPAHISON OF ANALYTICAL AND EXPERIMEMEAL RESULTS 
Raylelgh-Ritz Treatment in Four Degrees of Freedom 


When equation (l) for harmonic motion is treated in the Raylelgh- 
Ritz manner there exist in general four roots of flutter speed and fre- 
quency which satisfy the equations of equilibrium, although one or more 
roots may have no physical significance. For the configuration tested 
and analyzed, the critical (lowest) flutter speed root (designated root A) 
corresponds to en oscillation involving appreciable proportions of the 
rigid-body stability modes at a frequency which is a fraction of the wing 
first-bending natviral frequency. The next higher speed root (designated 
root B) corresponds to the conventional wing flexure-torsion flutter and 
only minute amplitudes of rigid-body motions are present. 

During the flight time after lavinchlng, several flutter pairameters 
were continually changing. Such varying parameters were; air density, 
weight and mass moment in_ pitch of the missile, euid the location of the 
missile center of gravity. These combined changes have only a moderate 
effect, an Increase of about 7 or 8 percent, on the calculated flutter 
speed (root A) from t = 0 to the time of wing failure. The calculated 
flutter speeds as functions of time are Included in figure 4. (The 
corresponding changes in flutter frequencies, amplitude ratios, and 
phasings are also small.) The intersections of calculated flutter speeds 
and flight speeds give the predicted flutter speeds, 520 feet per second 
for model D and 455 feet per second for model E, as shown in figure 4. 

The corresponding predicted flutter frequencies are 6.0 cycles per second 
for model D and 4.5 cycles per second for model E. Each configuration 
did not flutter, however, until its flight speed was 50 percent higher, 

670 feet per second for model D and 58 O feet per second for model E, as 
also shown in figure 4. Thus, the analysis was conservative in its pre- 
dictions of flutter speeds. 

With regard to the degree of conservatism of the analytically pre- 
dicted flutter speeds it should be pointed out that, even thbu^ the 
aspect ratio is 7 (including body intercept), the finite-span correction 
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to the flutter speed is large for small values of the reduced frequency- 
parameter k. For models D and E, the analytical value of k was approx- 
imately 0.027 which is small for flutter work. Reference 8 gives some 
results of flutter tests made in part to determine the effect of finite 
aspect ratio on the flutter of cantilever wings . The flutter was the 
hending-torsion type and occvirred mostly in a range of k from about 
0.15 to 0.30. According to the results of reference 8 for an aspect 
ratio of 6 to 8, the finite span causes an increase of rotighly 5 "to 10 per- 
cent in the flutter speed. For such a low -value of k as was obtained 
with models D and E, it is kno-wn that the finite span causes a consider- 
ably larger Increase In the flutter speed. Therefore, it appears that, 
if the proper finite-span corrections were made, an analytically predicted 
flutter speed wooild agree much better -with the experimen-tally determined 
one. Compressibility effects also Influence the results and can be taken 
into account. It is of Interest to note that an oscillation with k = 0.027 
has a period of II7 chords per cycle since, rt/k eq\ials the period. This 
■value of the period is well within the reuige of period treated in dynamic- 
stability work with rigid or undeformable aircraft. 

The amplitude ratios and phasing of the' -vario-us degrees of freedom 
associated with root A for model E when t = 5*2 seconds based on a wing- 
tip-bending amplitude of 1 inch are as follows; 


Degree of freedom 

Relative amplitude 

Relative phase angle, 
deg 

Wing-tip bending 

1.0000 inch 

0 

Wing-tip torsion 

.00456 radian 

179.1 

Rigid-body pitching 

.01679 radian 

186.1 

Rigid-body vertical 
translation 

.2457 inch 

186.0 


The parameters used were; mass moment of missile in pitch, 13.25 slug- 
feet^j weight of missile, 211 povinds; air density, 0.002260 slUg 
per cubic foot; and the parameters of the right wing of model E. The wlng- 
tip-torslon amplitude is one-quetrter of the body-pitching amplitude, and 
the vertical translation anplitude is one-quarter of the wing- tip-bending 
amplitude. Each degree of freedom is veiy nearly in phase or almost dia- 
metrically out of phase with other degrees of freedom, (on the basis that 
positive translation and wing bending are downward, whereas positive angu- 
lar displacement is leading edge or nose up) . Also of interest is the 
virtual identity of phasing of the two rigid-body modes, which fact, together 
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with the appropriate amplitude ratios, indicates that the missile is 
pitching effectively about a lateral sixls l 4.9 inches forward of the mis- 
sile center of gravity. All the analytical amplitude ratios and phase 
relationships agree with their experimental covinterparts (reported in 
section entitled "Experimental Results") within the experimental acciiracy. 

The anplitude ratios of wing flexure-torsion flutter (root B), based 
on a wing- tip-bending amplitude of 1 inch, axe minute for the two body 
degrees of freedom, whereas the wing- tip- torsion amplitude is O.O676 radian 
with a phasing of about -l8°. The ratio of body weight to wing weight 
ranged from 58 at time of launching to 51 at time of wing failure for 
model E; if this ratio were smaller, more relative body motion might be 
associated with root B. 


Effect of Various Binary and Ternary Combinations of Freedoms 

Analyses employing various binary and ternary combinations of degrees 
of freedom were made and the results were compared with those of the foxir- 
degree-of-freedom analysis. Results axe given in table II. The parameters 
used throughout are those of the right wing of model D and those listed at 
the top of table II. The four-degree- of-freedom analysis is designated 
type (a). 

Consider root A. The addition of the wing- torsion degree of freedom 
to an analysis not already including it corresponds to a reduction of 
torsional stiffness from an infinite value to the finite experimentally 
determined value. Even though the amplitude of wing-tip torsion in a 
fovir-degree-of- freedom analysis, type (a), is s mal l compared with body- 
pitching amplitude, addition of wing torsion to analysis types (c) or (e) 
so that they become type (d) or (a), respectively, causes a significant 
reduction of about 15 percent in the analytically determined flutter 
speed. These reductions indicate that it is potentially dangerous to 
exclude the wing torsion from an analysis of body- freedom flutter. The 
addition of the vertical- translation degree of freedom to analysis 
type (c) or (d), so that analysis type (e) or (a), respectively, is 
obtained, effects an increase of 50 percent in the flutter speed and 
demonstrates clearly that it is not realistic to exclude vertical trans- 
lation from an analysis treating body freedoms, at least for the present 
configuration. Thus, it seems clear that at least the four degrees of 
freedom employed in this report shovild be included in an analysis of sym- 
metrical body- freedom flutter. The wing-bending stiffness is of predomi- 
nant importance to the flutter speed for root A inasmuch as that flutter 
speed is virtxially directly proportional to the natiiral bending frequency, 
at least in the range of low ratio of bending frequency to torsion 
frequency . 

The flutter speed and frequency of root B, primarily wing bending- 
torsion flutter, are lanchanged regardless of which type of analysis, (a). 


NACA TN 5511 


15 


(b), or (d), is employed, and for root B tbe flutter speed is predominantly 
affected by the torsional stiffness in contrast with the predominant effect 
of bending stiffness on the flutter speed of root A. 

Effect of Relative Positions of Wing and Missile 
Center of Gravity 

Another variable whose effect was studied was the distance from the 
missile center of gravity to the wing leading edge. The flutter parameters 
used were again those of the right wing of model D and those listed in 
table II except for the fixed distance between the wing leading edge and 
missile center of gravity. The analysis was type (a); that is, it Included 
all four degrees of freedom. 

Figure 8 shows the two calculated roots A and B of flutter-speed 
coefficient as a function of the distance of the wing leading edge behind 
the center of gravity of the missile. The experimental flutter-speed 
coefficient for model D is inclvided in the figure at a wing position of 
5.38 semichords, even thovigh the experimental parameters differed slightly 
or moderately from those at the top of table II on which the ajialytical 
curves are based. If the paraineters of the ana lysis coincided with their 
experimental counterparts, the calculated curves would be affected by only 
a few percent. The in 5 >ortant point is that the experimental flutter speed 
is about 50 percent higher than the calculated one. It is pertinent to 
recall that the missile has no horizontal airfoil surface other than the 
wing and that the missile would become statically unstable with the wing 
very far forward near the center of gravity. Such a rigid-body static 
instability is not indicated in the figure. Two low-frequency failures 
of missiles which did have horizontal tails and which had the missile 
center of gravity at approximately one-eighth chord have been observed 
during rocket-propelled model tests made at the Langley Pilotless Aircraft 
Research Station at Wallops Island, Va. 

If analysis types (c) or (d) (see table II) which do not include the 
freedom of vertical translation were enq)loyed, it would be predicted, on 
the basis of the aerodynEunlc forces of two-dimensional inconq)ressible 
potential flow developed in reference 7^ that single-degree-of-freedom 
rigid-body pitching instability would be encountered over the approximate 
range of wing position shown in figure 8. This single-degree pitching 
instability was investigated analytically by Smllg in reference 9- The 

quantity la/jtpb^ of reference 9 (where Iq, represents the mass moment 
in pitch of the missile divided by the span) is about ten to eleven 
thousand for models D and E. The predicted borderline flutter speed 
(root A) , if no vertical translation were pemlltted, would be zero in the 
critical range of wing location, since no mechanical restraint of pitching 
exists; that is, the In-vacuo natural rigid-body pitching frequency is 
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zero. It appears, however, that such an unhappy possibility is precluded, 
even in two-dimensional incompressible flow, by the actual existing free- 
dom of vertical translation of the missile in flight, at least for the 
configuration treated. 

The two low-frequency failures of missiles observed during the rocket- 
propelled model tests mentioned previovisly will now be considered. Each 
of these models had a pair of unswept rectangular uniform model wings 
2k inches by 8 inches. For convenience, the models are designated models 6 
and 7- Tlie wings of model 6 were designed to be similBLr to, although 
subtly smaller than, the wings of models D and E of the present report, 
and the wing first-bending and first- torsion frequencies of model 6 were 
practically identical to those of model D. The leading edges of the 
wings of models 6 and 7 were located slightly ahead of the missile center 
of gravity, - 0.25 semi chord for model 6 and -0.38 semi chord for model 7* 

The complete models 6 and 7 hsui weights and mass moments in pitch approx- 
imately one- third as great as those of models D and E. Each model inclvided 
a horizontal tail having an area eqml to about 20 percent of the wing area, 
and a tail length of about one-quarter of the total tip- to- tip wing span. 

Thus, it may be seen that models 6 and T , except for the added hori- 
zontal tails, were dynamically fairly similar to a hypothetical model D 
or E which has had its wings moved forward near to the missile center of 
gravity. In the light of such observation, it is Interesting to examine 
figure 8 for a small negative value of the abscissa or wing position. 

In that region the speed of the flutter involving body motions is pre- 
dicted to be lower than the speed for wing bending- torsion flutter of a 
modified model D or E. 

All these clrc\amstances lea^ to the conjecture that the horizontal- 
tail volianes of models 6 and 7 were sufficiently large to prevent a 
static instability, but not sufficient, or the tails were ln 5 )roperly 
placed, to prevent a low-frequency body- freedom flutter. 


CONCLUSIONS 


The following conclusions may be drawn from flight tests and a 
mathematical analysis made to investigate a type of flutter Inyolving 
wing deformation and body motions of a particular configuration and 
certain modifications of that configuration: 

1. The possibility of a flutter type of dynamic instability 
involving appreciable proportions of rigid-body motions as well as wing 
structural deformations has been predicted analytically eind confirmed 
experimentally. The flutter was entirely different from conventional 
vlng bending-torsion flutter in that the frequency was about one-half 
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the wing firet-hendlng natiural frequency. The period was approximately 
100 chords per cycle which is well within the i^ange of period found in 
rigid-hody control- free stability work. 

2. A flutter einalysls based on two-dimensional incompressible flow 
provided a conservative (by 50 percent) prediction of the airspeed at 
which the flutter wovild occur for the configuration studied. 

5. Wing-bending stiffness rather than wing-torsional stiffness has 
the predominant effect on speed of the body- freedom flutter. 

4. At least four degrees of freedom, wing bending, wing torsion, 
body pitching, and body vertical, translation, should be included in an 
analysis of symmetrical flutter involving wing deformations and body 
motions . 


5. The analysis predicted that, if the wing were moved rearward from 
a position of coincidence of wing quarter chord and missile center of 
gravity (no horizontal tail present), the speed of flutter involving body 
motions wo\ild first decrease sharply xmtil a certain wing position rela- 
tive to the center- of-gravity location was reached and then would increase 
slowly as the wing was moved farther to the rear. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 28, 1950. 
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TABLE I 


WING PARAMETERS 


Parajueters 

Model D 

Model E 

Left wing 

Right wing 

Left wing 

Right wing 

Airfoil section 

NACA 65AOO6 

NACA 65AOO6 

NACA 65AOO6 

NACA 65AOO6 

Mcr 

0.84 

0.84 

0.84 

0.84 

1, in 

29.875 

29.875 

30 

30 



5 

5 

5 

5 

b, ft 

0.419 

0.419 

0.417 

0.417 

Xcg, percent chord . . . 

58.8 

37.6 

40.0 

40.94 

Xea^ percent chord . . . 

55.5 

55.2 

56.75 

57.5 

a 

-0.55 

-0.296 

-0.265 

-0.25 

a + Xa 

-0.224 

-0.248 

-0.20 

-0.1812 

1 

K (stnd.) 

24.6 

25.1 

51.11 

31.08 

r ^ 



0 .2255 

0.2105 

0.183 

0.1837 

fhi^ cps 

19.5 

20 

16 

16.5 

fh2, cps 

115.0 

115.5 

99 

106 

ft, cps 

154 

158.5 

118.5 

119 

fa, cps 

151 

157.5 

117.0 

117.5 

GJ, lb- in. 2 

3.56 X 105 

5.61 X 105 

2.552 X 105 

2.555 X 105 

El, lb-in. 2 

gv, 

2.76 X lo5 

0.05 

0.01 

2.95 X lo5 

2.51 X lo5 
0.04 

2.34 X lo5 
0.04 
0.01 

8a 

0.02 
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TABLE II 

MISSILE PABAMETERS 


[Missile center of gravity is 17 inches ahead of wing 
leading edge; missile mass moment in pitch is 
ll<-.8 slug-ft2; missile weight is 256 lb; air 
density p is 0.002272 slug per cu ft^ 


Type of 
analysis 

Degrees of 
freedom included 

(iL 

Flutter- speed 
coefficient, 
v/boJh 

Flutter 

frequency ratio, 
a>/o)h 

Root A 
(2) . 

Root B 
(5) _ 

Root A 

Root B 

a 

F 

T 

P 

V 

9.6 

22.5 

0.26 

5.40 

b 

F 

T 

- 

- 

— 

22.5 

— 

5.40 

c 

F 

- 

P 

- 

7.5 

— 

0.20 

— 

d 

F 

T 

P 

- 

6.k 

22.5 

0.17 

5.40 

e 

F 

- 

P 

V 

11.2 

— 

0.287 

— 


^Designation of degrees of freedom 
F - wing flexure (first bending) 

T - wing torsion (first torsion) 

P - rigid-body pitching 
V - rigid-body vertical translation 

^oot A - critical body-freedom root 
5Root B - critical wing bending- torsion root 


Sta. 

0 58.7 76.0 92.7 



: NACA 65AOO6 

: NACA 65-009 


Figure 1 .- Layout of flutter test vehicle, model E. All dimensions are 

in inches. 
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Moment of inertia I about center of gravity, slug-ft 
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Figure 2.- Variations of model E parameters with time 
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Figure 3.- Model E in launching position. 
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Figure 4 .- Velocity plotted against time. 
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density, p, slugs/cu ft 
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Figure 5-- Air density and Mach number variations with time. 
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Figure 6.- Variations of longitudinal acceleration with time. 



Figure 7.- Portion of telemeter record of model E 
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Distance of wing leading edge behind model center of gravity in wing semichords 


Figure 8.- Variation of flutter-speed coefficient with wing position. 
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